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Abstract 

Starting from constrained Skyrme- mean- field calculations, the low-energy excita- 
tion spectrum of 16 O is calculated by configuration mixing of particle- number and 
angular-momentum projected mean-field states in the framework of the Generator 
Coordinate Method. Without any adjustable parameters, this approach gives a very 
good description of those states and their transition moments that can be described 
with our restriction to axially and reflection-symmetric shapes. The structure of 
low-lying + states is analyzed in terms of self-consistent Op-Oh, 2p-2h, and 4p-4h 
Hartree-Fock states. 
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1 Introduction 

The doubly-magic oxygen isotope 16 has always attracted much interest in 
nuclear structure physics. Its first excited + state has been interpreted by 
Morinaga [1] as an evidence for shape coexistence. Multi-nucleon transfer re- 
actions have lead to the conclusion that the presence of a state corresponding 
to a deformed shape is due to excitations of both proton and neutron pairs 
across the N and Z = 8 closed shells. Several other excited states at low ex- 
citation energies can be interpreted as due to multi-particle-hole excitations 
(see the discussion in the Wood et al. review paper [2]). The first excited + 
state is the head of a rotational band. 

Shell model and mean-field calculations were already performed in the sixties 
[3,4,5] with a description of this + state based on a deformed 4p-4h configu- 
ration. The cluster model [6] , using an empirical a-a interaction, qualitatively 
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explains the deformed structure of the excited states. To go beyond a quali- 
tative understanding of the 16 spectrum appeared however very quickly as 
a difficult challenge. Mean-field calculations with more modern effective in- 
teractions [7,8,9] clearly showed that correlations beyond the Hartree-Fock 
(HF) approach have to be included to describe successfully the first excited 
+ state. Two kinds of correlations were mainly invoked: pairing correlations, 
and restoration of rotational symmetry which both leave the spherical config- 
uration unchanged and lower the energy of the deformed configurations. 

More recently, shell model calculations mixing (0+2+4) fiw excitations [10,11] 
have confirmed the schematic model of Brown and Green [4] , and showed that 
the large quadrupole moment of the deformed configuration appears naturally 
in a full shell model diagonalization. It has also been possible in such calcu- 
lations to test how well 16 O can be considered as an inert core as it is usually 
done in shell model calculations. The Op-Oh component of the ground-state 
wave function was found to have a weight of the order of 40 to 50% only, with 
nearly equal importance of the 2p-2h components. 

The question of a correct description of the 16 O spectrum starting from a 
mean-field approach is thus still an open question. This nucleus is an an- 
chor point where many approaches to the nuclear many-body problem can be 
compared. Indeed, besides the already mentioned mean-field and shell model 
calculations using effective interactions, there are also calculations based on 
so-called "realistic" nucleon-nucleon interactions (introduced as effective po- 
tentials to describe the phase shifts in nucleon-nucleon scattering), either al- 
ready performed for 16 [12], or to be expected for the near future [13]. 

In this paper, we present an application of the method introduced by Valor 
et al. [14], which performs a configuration mixing of projected self-consistent 
mean-field states, to calculate the low-energy spectrum of 16 0, and to analyze 
its structure in terms of np-nh Hartree-Fock states. 



2 The Method 



The starting point of the method is a set of BCS states \q) generated by con- 
strained Skyrme-mean-field calculations. We have used the Skyrme parametriza- 
tion SLy4 [15] together with a like-particle T = 1 density-dependent zero-range 
pairing interaction [16]. The pairing strength adjusted in [17] has been reduced 
from V = -1250 MeV fm 3 to V = -1000 MeV fm 3 for both protons and neu- 
trons, as in [14]. The pairing active space is limited by a soft cutoff at 5 MeV 
above and below the Fermi energy [18]. This combination of effective inter- 
actions has been proven to be very successful in the description of a large 
number of experimental data all over the chart of nuclei. 
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The set of mean-field wave functions \q) is generated by calculations with a 
constraint on the axial quadrupole moment Qq. To avoid a collapse of pair- 
ing correlations along the constraining path, the approximate variation-after- 
projection Lipkin-Nogami (LN) method is used. The Lipkin-Nogami prescrip- 
tion to evaluate an energy correction due to particle number projection is 
known to have deficiencies when pairing correlations are weak [19]. We do not 
make use of this correction, however, but introduce the LN method solely to 
ensure that pairing correlations, although weak, are present even in the spher- 
ical configuration. With that, the LN method provides BCS wave functions 
which are a fair starting point for an exact projection on particle number. 
Since these projected BCS wave functions are used further as basis functions 
for a configuration mixing, the deficiencies of the LN method should not affect 
the results. 

The most important symmetries broken by the mean-field approach are re- 
stored after variation by standard projection techniques [14]. For each value 
of the quadrupole moment, the wave functions are projected simultaneously 
on angular momentum and on neutron and proton particle numbers, decom- 
posing a given intrinsic BCS state \q) into wave functions | Jq) corresponding 
to several values of the angular momentum. For each value of J separately, 
these sets of non-orthogonal states \Jq) are then mixed as a function of the 
quadrupole moment to give the final wave functions \Jk), which correspond to 
collective states of the nucleus. A discrete number of values of the quadrupole 
moment are considered, but in such a way that the results do not depend on 
the discretization and are equivalent to a mixing on a continuous variable as 
in the generator coordinate method (GCM) [20,21]. Among the configuration- 
mixed states \Jk) so obtained, the physically interesting ones are the ground 
state (J = 0, k — 0) and the few excited states which can be described 
by an axial quadrupole collective mode. States corresponding to modes that 
require a breaking of reflection symmetry (like a- 12 C configurations or the 
tetrahedron-like 4-a configuration) are not included in our description of 16 O. 



3 Projection 

The energy curves are shown in figure 1. The mean-field curve, projected on 
particle numbers N = Z = 8 only, shows, as expected, a deep minimum for 
the spherical configuration and no secondary minima. Two inflection points, 
however, are present at large deformations (80 and 150 fm 2 ). Looking at the 
single-particle level scheme plotted for neutrons in figure 2, one sees that 
these points correspond to level crossings. At a /?2 value around 0.8 (or a 
mass quadrupole moment equal to 80 fm 2 ), a rf 5 / 2 level crosses the pi/2 level 
and is located below the Fermi level for larger quadrupole moments. Since 
due to time-reversal invariance, the levels are doubly degenerate and since 
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Fig. 1. Particle-number projected ("mean field") and particle-number and angu- 
lar-momentum projected potential energy curves for 16 0. All curves are drawn ver- 
sus the (intrinsic) deformation of the BCS state which is projected. The upper axis 
gives the mass quadrupole moment Qo = a/167t/5 (r 2 Y 2 o), the lower the dimen- 
sionless deformation parameter j3 2 = ^ {r 2 Y 20 ) / {2>R 2 A) with R = 1.2 A 1 / 3 fm. The 
various shapes along the paths are indicated by the contours of the total density at 
po = 0.07 fm~ 3 . 

neutrons and protons have a very similar behavior in this N = Z nucleus, 
the mean-field wave functions for j3 2 values larger than 0.8 correspond to a 
deformed 4p-4h configuration. Owing to pairing correlations, however, all the 
levels around the Fermi energy have sizable occupation probabilities. In the 
same way, for a f3 2 value around 1.8, a crossing occurs between a state coming 
down from the pf shell and a p state. Note that for the protons the level cross- 
ings are at slightly smaller deformation than for the neutrons, a consequence 
of the Coulomb interaction. For extremely deformed shapes, above f3 2 = 4.0, 
the mean-field configuration corresponds to a deformed 8p-8h configuration 
which is associated to a chain of four a particles. 

The energy curves projected on angular momentum, from J = to J = 6, are 
also shown on figure 1. While the spherical mean- field state has a J = com- 
ponent only, the deformed configurations can be projected on several angular 
momenta. Their J = components have always energies lower than the mean- 
field energy. Since for each quadrupole moment the mean-field energy is the 
weighted sum of the projected energies, figure 1 shows also that components 
with J > 6 are dominant beyond f3 2 ~ 2.5. 



4 



Qo (fm 2 ) 

100 200 300 400 500 




ft 

Fig. 2. Single-particle spectra for neutrons in 16 0. Solid (dotted) lines represent 
levels with positive (negative) parity. 

The J = curve is extremely flat up to the first single-particle level crossing. 
Beyond this point, the gain in energy due to angular momentum projection 
is approximately constant and of the order of 7 MeV. For angular momenta 
equal to 2 to 6, a well-developed minimum is obtained for mean-field states 
with an axial quadrupole moment around 75 fm 2 which correspond to the 
deformation at which single-particle level crossing occurs. The shape of the 
mean-field density is also changing abruptly at this point, as can be seen on 
figure 1 and is reflected in the hexadecapole deformation. For a quadrupole 
moment value around 100 fm 2 , the fa parameter, as defined in [22], changes 
sign. For still more deformed intrinsic wave functions, the projected energies 
vary very slowly, with soft minima obtained for fa values around 1.5. 



4 Configuration Mixing 

The configuration mixing of mean-field wave functions as performed in this 
work has several goals. The particle-number projection removes unwanted con- 
tributions coming from states with different particle numbers. The angular 
momentum projection separates the contribution from different angular mo- 
menta and generates wave-functions in the laboratory frame which provide 
transition probabilities and spectroscopic moments without further approxi- 
mations. Finally, the variational configuration mixing with respect to a col- 
lective coordinate, the axial quadrupole moment in this work, removes the 
contributions to the ground state coming from collective vibrations and pro- 
vides the excitation spectrum corresponding to this mode. 
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State E'expt ^calc ^rms Qs Qo{s) (h{s) 

(MeV) (MeV) (fm) (e fm 2 ) (e fm 2 ) 



07 


0.00 


0.00 


2.75 










6.05 


6.03 


2.90 










6.92 


7.09 


2.91 


-11.7 


41.0 


0.74 




12.05 


12.45 


3.32 








2+ 


13.02 


13.18 


3.22 


-22.9 


80.2 


1.4 


o 4 + 




20.73 


3.85 








n 




23.06 


3.82 


-44.4 


155 


2.8 



Table 1 

Experimental (-E e xpt) and calculated (Scale) excitation energies for the low-lying + 
and 2 + states in 16 0. Only states that can be described by our calculation are indi- 
cated. Also given are the calculated proton rms radii r rms , spectroscopic quadrupole 
moments Q s (in the laboratory frame), proton quadrupole moments Qo(s) and /^(s) 
deformation parameters in the intrinsic frame. Experimental excitation energies are 
taken from [23]. 

Results for the lowest states obtained in the configuration mixing calculation 
are given in table 1 and are compared with the experimental excitation en- 
ergies. The calculated excitation spectrum is compared to the experimental 
data in figure 3 and to the mean-field energy curve in figure 4. Only pro- 
late configurations have been included. The projected oblate configurations 
have an overlap close to 1.0 with the spherical configuration which makes the 
configuration mixing calculation numerically unstable. The rms radii and de- 
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Fig. 3. Comparison between the experimental and theoretical low-lying + and 2 + 
states in 16 O. 
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Fig. 4. Spectrum of the lowest + (circles) and 2 + (triangles) states, plotted at a 
deformation corresponding to the average deformation of the intrinsic wave func- 
tions. These averages are determined using the weight of each mean- field state in 
the collective wave functions. Also shown are the angular momentum projected en- 
ergy curves and the (particle-number projected only) "mean-field" curve (thin black 
line) . 

formations given in table 1 are calculated from the point-proton distribution 
and do not include the corrections that have to be introduced to calculate 
charge radii that can be compared with the experimental data, see e.g. [24]. 

With our method, one obtains automatically the spectroscopic quadrupole 
moments Q s in the laboratory frame. The charge quadrupole moment in the in- 
trinsic frame Qo(s) is then calculated with the relation Qo(s) = —Q s (2 J + 3)/ J 
assuming K = states [25]. The dimensionless deformation parameter fa is 
related to Q (s) by fa(s) = ^5/(16tt) x AttQ (s) /(3R 2 Ze 2 ). 

The assignment of the calculated and 2f states to experiment in table 1 is 
tentative. Note that the restriction of our calculation to reflection-symmetric 
configurations does not permit to describe states with negative parity or whose 
intrinsic structure is supposed to be asymmetric. For instance, the ct- 12 C con- 
figuration assigned to the experimental 2 + state at 9.8440 MeV [23] cannot 
be represented by our model. For the other states, the agreement between the 
experimental and theoretical energies is excellent. 

The correlations introduced by the configuration mixing decrease the en- 
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Fig. 5. GCM wave functions \Jk) = J2 q 9k(<l)\Jq} of the four lowest + states. 

ergy of the 0+ state with respect to the HF ground state by 2.3 MeV. This 
gain in energy, although significant, is lower than in several of our previous 
configuration-mixing calculations. In 208 Pb, for instance, removing quadrupole 
vibrations from the mean-field ground state by configuration mixing without 
angular momentum projection brings already 2.0 MeV of extra energy [26]. 
In 24 Mg, restoration of symmetries and the inclusion of the contribution from 
quadrupole vibrations decrease the mean-field ground-state energy by more 
than 5 MeV [14]. 

The mixing coefficients Qk{q) for the four lowest + states are plotted in figure 
5. The ground-state wave function is spread over a large range of deformations, 
with nearly equal weight up to a j3 2 ~ 0.8. Its rms radius r rms = 2.75 fm is 
therefore significantly larger than the 2.68 fm value of the spherical HF state. 
The maximum weight for the first excited state is obtained around fa equal to 
1.0, but with significant components at lower quadrupole moments. The third 
and fourth + states are shifted to configurations with larger deformations. 
As can be seen from the mean deformation of their mean-field components 
(see Fig. 4), and the shapes displayed in Fig. 1, none of them corresponds 
to a chain of 4 a particles which occurs at even larger deformation around 
fa « 4.0. 

The B(E2) values for transitions from the 2\ to the ground state and to the 
first excited + state are well described by our model. The error is larger for 
the transition from the l\ state to the ground state. This is not surprising 
since important configurations like a— 12 C which are not included in our model 
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transition 


B(E2) T /(e 2 fm 4 ) 


Qo(t) 


lh(t) 


transition 


M(M))/fm 2 


expt. 


calc. 


(e fm 2 ) 


expt. calc. 


2+^0+ 


42 ± 1.4 


38 


19.7 


0.36 


0+-0+ 


3.55 ±0.21 5.735 


2+-0+ 


21 ±7 


4.7 


6.9 


0.12 


0+^0+ 


4.03 ± 0.09 0.690 


2+^0+ 


370 ±4 


241 


49.2 


0.89 






2+^0+ 


42 ±7 


74 


27.3 


0.49 






2+-0+ 




590 


77.0 


1.39 







Table 2 

Reduced transition probabilities B{E2) |, calculated transition quadrupole mo- 
ments Qo(t) = y/l§itB(E2) f /5, corresponding deformation parameters fh(t), and 
monopole matrix transition elements M(E0) = (f\ ^ p fp\i) between low-lying states 
in 16 0. The experimental data are taken from [23]. 

should play a role at the excitation energy of the 2^ state. For a similar reason, 
the monopole transition probabilities are in better agreement with experiment 
for the first excited + state than for the second. 

Table 2 lists the experimental value from [23]. The values adopted by Raman 
et al. [27] for the B{E2) |= (40.6 ± 3.8) e 2 fm 4 of the 2f -> Of transition is 
even in better agreement with our result. The Qo(t) of the transitions within 
a rotational band, i.e. 2f — > 0^ and 2\ — > O3 , are close to the intrinsic 
quadrupole moments Qo(s) of the 2 + states involved as given in table 1. 



5 Analysis of the Collective States 

Analyses of the shell model wave functions obtained in a (0+2+4) fuv active 
space [10,11] have shown that the 16 ground-state wave function is by far 
not a pure Ohu state and that the breaking of the closed shells at iV and Z 
equal to 8 is large. The Ohu component of the shell-model wave function has a 
weight of the order of 50 % only, with components of comparable importance 
in the 2 fko subspace. The first excited + state has an energy very close to 
the first experimental + and is composed to 90% of a 4p-4h configuration. 

It is not evident to perform a similar analysis for self-consistent wave functions. 
In our model, there is no spherical oscillator basis on which the individual wave 
functions are expanded. However, it is important to find a common language 
with shell-model calculations and to interpret a deformed mean-field basis in 
terms of a spherical basis. In the framework of self-consistent mean-field mod- 
els, a natural choice for the reference state is provided by the self-consistent 
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Hartree-Fock (HF) Op-Oh ground state, which is of course spherical for the 
closed-shell nucleus 16 0. To build np-nh excitations in an HF framework re- 
quires some care. They could be constructed non-self-consistently on top of the 
spherical HF state by just changing the occupation of the single-particle levels 
below and above the Fermi surface. However, such excitations are not unique 
and the resulting wave function is, in general, not spherical. A much more nat- 
ural procedure is to consider the deformation also in the np-nh HF states and 
to construct them self-consistently by solving the HF equations with appropri- 
ate occupation numbers of the single-particle states. These 4p-4h and 8p-8h 
states have deformations corresponding to the deformed shell closures that can 
be seen around fa ~ 1-4 and /?2 ~ 4.5 in figure 2. Each fully self-consistent np- 
nh HF state | np-nh) defines a basis of deformed single-particle states different 
from the other excitations. As a consequence, the various |np-nh) are not or- 
thonormal. However, since of their very different structure (deformation and 
np-nh excitations), their overlap is extremely small and can be neglected. 

Since they are deformed, the np-nh HF states are not eigenstates of angular 
momentum and projected np-nh states | Jnp-nh) are constructed by symmetry 
restoration. 

We analyze the np-nh content of the collective states \Jk) by calculating the 
amplitude | (Jk\ Jnp-nh) | 2 . Note that our procedure is unambiguous and can 
be used to analyze any state constructed from mean-field wave functions. It 
permits to quantify the differences between simple mean-field configurations 
and the fully projected configuration mixing wave functions. The way np- 
nh excitations are constructed is not the same as in shell model calculations 
and the single-particle wave functions do not have the same analytical form 
(the self-consistent wave functions discretized on a mesh have a much better 
asymptotic behavior than oscillator wave functions). Nevertheless, our proce- 
dure will give us some insights on the relation between our wave functions and 
shell model ones. 

The properties of the HF and projected HF states are given in table 3 where 
we also indicated the weight of the projected HF states in each of the four 
lowest collective + states. 

From figure 6, one clearly sees that there is a correlation between the de- 
formation of the projected mean-field state and its overlap with the np-nh 
configurations. One sees also from table 3 that the deformation of the np-nh 
HF states increases rapidly with the number n of excitations. 

With about 75 % weight of the Op-Oh HF state, see table 3, the breaking of the 
closed shell character of the Of ground-state wave function is not as large as 
in shell model calculations, but still sizable. Note that a part of the admixture 
of higher shells obtained in shell-model calculations is due to the fact that 
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state E HF E mF Q (3 2 \(Jk\ Jnp-nh)\ 2 





(MeV) 


(MeV) 


(fm 2 ) 




ot 


o 2 + 


Os 


O4 


Op-Oh 


0.0 


0.0 





0.00 


0.737 


0.195 


0.006 


0.0001 


2p-2h (p) 


14.92 


11.96 


41 


0.37 


0.094 


0.098 


0.018 


0.003 


2p-2h (n) 


15.20 


12.09 


41 


0.37 


0.099 


0.098 


0.019 


0.002 


4p-4h 


15.95 


11.26 


119 


1.08 


0.027 


0.212 


0.193 


0.003 


8p-8h 


38.03 


31.48 


466 


4.21 


0.000 


0.000 


0.000 


0.034 



Table 3 



Excitation energy of the HF (Ehf) an d J = projected HF states (Ejbf) with 
respect to the (Op-Oh) HF ground state. The mass quadrupole moment Qq and 
the deformation f3 2 of the self-consistent Op-Oh, 2p-2h, 4p-4h, and 8p-8h HF states 
are also given. The last four columns give the weights \(Jk\Jnp-nh)\ 2 of their J = 
component in the four lowest + states | Jk) obtained from the configuration mixing. 

the "real" single-particle potential is not an harmonic oscillator as assumed 
in a shell model basis. Therefore it is not surprising that we find a smaller 
admixture of higher shells as we start with single-particle wave functions with 
a more realistic asymptotic behavior. 

For the 4p-4h configuration, where a pair of protons and neutrons is excited 
from the lpi/ 2 to the 2d 5 / 2 levels, a minimum in the HF energy is obtained for 
a mass quadrupole moment of Q — 119 fm 2 , or (3 2 = 1-08. This 4p-4h config- 
uration corresponds to the deformed N = Z = 8 shell closure visible at about 
the same deformation in figure 2. The excitation energy of this deformed HF 
state (with respect to the spherical Op-Oh HF state), is equal to 15.95 MeV, 
far above the experimental value of 6.03 MeV for the first excited + state. 
The restoration of rotational symmetry by projection on angular momentum 
J = decreases the calculated excitation energy to 11.25 MeV, still too high 
compared to experiment. To complete the analysis, we have also constructed 
2p-2h and 8p-8h configurations. The only two 2p-2h states that we have con- 
sidered are obtained by promoting either a pair of protons or a pair of neutrons 
from the pi/ 2 shell to the d 5 / 2 shell. As we include the Coulomb interaction 
self-consistently, these two states have slightly different energies. The 8p-8h 
configuration corresponds to the occupation of the sd and pf states which are 
below the Fermi level at very large deformations (see figure 2). 

The Op-Oh HF state is the dominant component of the collective ground state. 
If one adds the weights of the Op-Oh, 2p-2h and 4p-4h states, one obtains a 
nearly complete description of the ground state, the missing components being 
probably 2p-2h states not included in our analysis. One should, however, take 
into account that the self-consistent np-nh states are not orthogonal and have 
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Fig. 6. Weight \(Jq\ Jnp-nh)\ 2 of the J = projected spherical Op-Oh and deformed 
2p-2h, 4p-4h, and 8p-8h HF states | Jnp-nh) in the J = projected paired mean- field 
states \Jq). The weight of the two different 2p-2h states is very similar, therefore 
just one is drawn. 

small, but non-vanishing overlaps. With less than 3%, the contribution from 
the pure HF 4p-4h configuration to the ground state is small and lower than 
the shell-model values (between 4 and 12% depending on the interactions). 

The excited state has a much more complicate structure, without any 
dominant configuration. As can be seen in table 3, the weights of the Op-Oh, 2p- 
2h and 4p-4h configurations are quite close. These four components, however, 
represent only about 60% of the full state. Looking to figure 5, one can see 
that the spreading of this state extends up to deformations of the order of 200 
fm 2 , where the overlap between the np-nh HF states and the projected BCS 
states are small (see figure 6). Many other np-nh configurations are therefore 
necessary to obtain a full description of this 0^ state. 

None of the collective states is dominated by the deformed N = Z = 8 shell 
closure. The excited O3 state has also a small 4p-4h component and is con- 
structed on mean-field states with larger deformation than the 4p-4h HF state. 

As can be seen from table 3, the 8p-8h state plays no role for the low-lying 
+ states. The first state in which it has a large weight of 0.48 is the 0^j~ 
state. The density distribution of the 8p-8h state is close to that of a four-a- 
chain configuration. Its 38.0 MeV excitation energy is lowered to 32.2 MeV 
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by projection on angular momentum and to 32.0 MeV by the configuration 
mixing. It is thus very large. This is, however, not inconsistent with the fact 
that there is up to now no convincing evidence for the presence of such a state 
in the low-energy spectrum of 16 O. 



6 Summary and Conclusions 

This study of the doubly-magic nucleus 16 has demonstrated the descrip- 
tive and predictive power of methods based on a self-consistent mean-field 
approach including correlations. In particular, the gain in energy of the pure 
4p-4h HF state, that is brought first by projection on angular momentum, 
and then by mixing on the axial quadrupole moment, is impressive. This pro- 
cedure does not only permit to obtain the first excited + state at the right 
energy, but also other excitations and transition matrix elements which are in 
reasonable agreement with available data. Somewhat surprisingly, the 4p-4h 
HF state is evenly spread over the first and second excited + states, leaving 
none of them associated with a simple 4p-4h configuration. 

Our results are not very sensitive to the parametrizations of the mean-field and 
pairing interactions that are used. The Skyrme interaction SLy4 used here was 
adjusted at the mean-field level to binding energies, charge radii, and nuclear 
matter properties of a few nuclei, among them 16 0. The correlations taken 
into account in our study increase the binding energy by 2.3 MeV compared 
to the spherical HF state, which leads now to a slight over-estimation of the 
16 O ground state-energy. Similarly, the proton rms radius is slightly increased. 
Although it is gratifying to see that the change in total energy is below 2% 
while the change in radii stays below 3%, it is clear that a readjustment of the 
effective interaction will have to be done when a larger experience will have 
been obtained on the effect of correlations. 

We have included the most important symmetry restorations to permit to 
compare our results directly with experimental data in the laboratory frame 
of reference. Our calculations demonstrate that this is a key to a successful 
quantitative description of the low-energy states when starting from a mean- 
field approach. Some ingredients, however, are still missing in our model. There 
are still symmetries broken by our approach. Our wave functions break trans- 
lational and Galilean invariances. The implicit assumption is made that these 
broken symmetries bring similar errors on all collective wave functions. To 
check this hypothesis in our model is unfortunately still beyond numerical 
possibilities. Projection on isospin seems to be the next natural step to enlarge 
the predictive power of the method even further. We have also not considered 
proton-neutron pairing correlations, although 16 is an N = Z nucleus. Devel- 
opments in the understanding of these correlations and their modeling are still 
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necessary. Concerning shape degrees of freedom, reflection-asymmetric config- 
urations like 12 C+ 4 He still have to be incorporated to describe several states 
of astrophysical interest. We have also restricted ourselves so far to axially 
symmetric shapes. Including triaxiality might alter the results for high-lying 
excitations which spread differently into the f3-^ plane. 
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